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Stark-induced amplitudes for the 6P1/2 — 7Pi/2 transition in Tl I are calculated using the rela- 
tivistic SD approximation in which single and double excitations of Dirac-Hartree-Fock levels are 
summed to all orders in perturbation theory. Our SD values Q5 = 368 Oq and |/3s|= 298 ag ^-re 
in good agreement with the measurements qs=377(8) Oq and Ps = 313(8) Oq by D. DeMille, D. 
Budker, and E. D. Commins [Phys. Rev. A50, 4657 (1994)]. Calculations of the Stark shifts in the 
6P1/2 — 7P1/2 and 6P1/2 — 7Si/2 transitions are also carried out. The Stark shifts predicted by our 
calculations agree with the most accurate measured values within the experimental uncertainties for 
both transitions. 
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PACS numbers: 31.15.Ar, 31.15.Md, 32.10.Fn, 32.70.Cs 



I. INTRODUCTION 

The 293 nm 6P1/2 — IP1/2 transition in thallium 
has been studied extensively both experimentally and 
theoretically because of its connection to atomic par- 
ity nonconservation (PNC). Although 6P1/2 — TPi/2 is 
nominally a magnetic dipole (Ml) transition, there is 
also an electric dipole (El) component arising from 
the weak interaction mediated by Zq exchange be- 
tween the nucleus and bound electrons. The weak 
El component of the 6P1/2 — 7Pi/2 transition is cal- 



ibrated using the Stark-induced amplitude j3s- Mea- 
sured values of both 0:5 = 2.01(4) xlO~^/iB cm/V and 
Ps = 1.67(4) X 10- Vb cm/V for the 6P1/2 - 7Pi/2 tran- 
sition in thallium were reported by DeMille et al. [ij 
and found to be in substantial disagreement with ear- 
lier measurements as — 1.31(6) xlQ-^/is cm/V and 
Be = 1.09(5) xlO-^/z_B cm/V by Tanner and Commins 
01 ■ It was pointed out in Ref. [3 that both the Stark 
shift and the Stark-induced amplitude could be measured 
in the same transition in which PNC is measured y. 

High-precision measurements of the Stark shift within 
the 378 nm 6Pi/2~7Sii2 El transition in atomic thallium 
were recently reported by Doret et al. |3]- The result 
Avs = -103.23(39) kHz/(kV/cm)2 had higher accuracy 
by factor of 15 than earlier measurements [l|, |a]. The 
earlier value for the 6P1/2 — 751/2 transition from Ref. |ll| 
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was Ai/s = -112(6) kHz/(kV/cni)2. 

In the present paper, we carry out relativistic all-order 
calculations of Stark-induced amplitudes as and /3s for 
the 6P1/2 — 7Pi/2 transition as well as Stark shifts within 



both 6P1/2— 75i/2 and 6Pi/2~7Pi/2 transitions in atomic 
thallium. The calculations are carried out using the rel- 
ativistic SD all-order method in which single and double 
excitations of Dirac-Hartree-Fock (DHF) wave functions 
are summed to all orders in perturbation theory. Re- 
cently, lifetimes. El, E2, and Ml transition rates, hy- 
perfine constants, and excitation energies of the nPj, 
nSi/2, and nDj states in neutral thallium were evalu- 
ated by Safronova et al. 6| , using both relativistic many- 
body perturbation theory (MBPT) and the relativistic 
SD method. The SD calculations were found to be in 
excellent agreement with the best available experimental 
data. 



II. STARK-INDUCED AMPLITUDE OF THE 

6P1/2 - 7Pi/2 TRANSITION IN THALLIUM 

Stark-induced scalar and vector polarizabilities as and 
(3s for transitions in Na, K, Rb, Cs, and Fr were calcu- 
lated in the SD approximation by Safronova et al. Q- 
Following the procedure used in [3|, we treat Tl I as a 
one-valence-electron atom with a Hg-like core and eval- 
uate the Stark-induced scalar and vector polarizabilities 
as and f3s following the procedures used in Q for alkali- 
metal atoms. 

The scalar and vector polarizabilities as and Ps for 
the transitions between the ground state 6P1/2 and the 



TABLE I: The contributions to Stark-induced scalar polarizability as for the 6pi/2 — 7pi/2 transition in Tl I. The corresponding 
energy differences and electric-dipole reduced matrix elements are also listed. All values are given in a.u. 



Contribution 


nlj 


^(>Pl/2 ~ ^-n-lj 


E^Pl/2 - 'E'nij 


•^7Pi/2,nij 


■^n!j.6Pi/2 


as 


aT""{nS) 


75i/2 


-0.120640 


0.035004 


6.016 


-1.827 


-37.2 




8<5'i/2 


-0.176539 


-0.020895 


-6.215 


-0.535 


-29.7 




9Si/2 


-0.196680 


-0.041036 


-1.274 


-0.298 


-1.9 




lOSi/2 


-0.206387 


-0.050743 


-0.656 


-0.200 


-0.5 


arinS) 












-2.0 


«r"(ni33/2) 


6D3/2 


-0.164565 


-0.008921 


-10.703 


-2.334 


492.0 




7O3/2 


-0.191418 


-0.035774 


4.821 


-1.101 


-29.4 




8^^3/2 


-0.203543 


-0.047899 


2.377 


-0.672 


-6.9 




903/2 


-0.210040 


-0.054396 


1.535 


-0.476 


-2.8 


ar'(nD3/2) 












-14.1 


SaT' 












0.4 


Total 












368 



excited state 7Pi/2 in Tl I are calculated as 

as = ^[la{nSi/2) - Ia{nD:i/2)\ , (1) 



/3s -E 



Ip{nSi/2) + -If3{nD:i/2) 



where 



Ia{nlj) 

__ 1 
~ 6 



(2) 



EePl/2 - ^nlj -E^7Pi/2 ~ ^nlj 



and 



_ 1 
~ 6 



Z7Pi/2,nljZnljfiPi/2 



^7Pi/2,«ii^nii,6Pi/2 



E7P,, 



E, 



nlj 



EqP 



E.^ 



■nlj 



The quantities Z^v in the above equations are electric- 
dipole reduced matrix elements. The aU-order SD matrix 
elements are calculated as Q 



^IDV 



Zi-) 



zw 



^{l + K^){l + N,) ' 



(3) 



where Zwv is the lowest-order (DHF) matrix element and 
the terms Z^^\ k — a- ■ -t are linear or quadratic func- 
tion of the SD excitation coefficients. The normalization 
terms iVu, are quadratic functions of the excitation co- 
efficients. The SD matrix elements include all MBPT 
corrections through third order together with important 
classes of forth- and higher-order corrections. 

The calculation of the as is divided into three parts: 



as 



„ main 



I „ tail I r„ core 



(4) 



where a™^" is the dominant contribution from states 
near the valence state, (5a|?''° is the contribution from 



core-excited autoionizing states, and ag is the remain- 
der from highly excited one-electron states. Thus, we 
write. 



10 



„ main 
"5 



^Ia{nSi/2)'-^Ia{nD^/2), (5) 

n—7 n=& 

6 5 



n=l 

N 



71=3 

N 



ataii = E -^"('^'^1/2)- E -^"("^3/2), 

n=ll n=10 

where N is the number of the finite basis set states. The 
calculation of /3s is conducted in the same way. 

We use B-splines [^ to generate a complete set of DHF 
basis orbitals for use in the evaluation of all electric- 
dipole matrix elements. Here, we use here A^ = 50 
splines (compared with 40 used in [3) for each angu- 
lar momentum. The basis orbitals are constrained to 
a cavity of radius i? = 90 (a.u.). The size of the cav- 
ity is taken to be large enough to fit all of the states 
needed for the calculation of the main terms for all of 
the polarizabilities calculated in this work. Furthermore, 
we use the Breit-Dirac-Hartree-Fock (BDHF) approxi- 
mation here instead of the Dirac-Hartree-Fock (DHF) 
approximation used in Ref. 6] because we found that 
the Breit contributes about 1% to the values of as and 
13s at the DHF level. The BDHF approximation includes 
the the one-body part of the Breit interaction in the DHF 
equation (for more detail see Refs. [^illfl)- The one-body 
part of the Breit interaction is also included in the equa- 
tions that generate DHF basis orbitals. In fact, all of 
the calculations in this work, including the calculations 
of the tail contributions, are done with same basis set. 

The sums over n in Eqs. (^12) converge rapidly; there- 
fore, only main term contributions have to be calculated 
accurately. We calculate main terms for as and (3^ using 
SD matrix elements and experimental energies 11]. The 



TABLE II: The contributions to Stark-induced vector polar- TABLE III: The SD all-order Stark-induced scalar and vec- 
tor polarizabilities for the 6P1/2 — 7Pi/2 transition in Tl are 
compared with measurements by Tanner and Commins Q - 
a and DeMille et al. ^l| - b. 



izability /3s(a.u.) for the 6P1/2 


— 7P1/2 transition in Tl I. 


Contribution 


nlj Ps 



/3J?"'"(nS) 



P's'^'^nS) 



PT"^{nDs/2) 



ops 
Total 



751/2 

8<S'i/2 
95i/2 

10^1/2 



6D; 

7D: 



3/2 

'3/2 



8-D3/2 
91)3/2 



-67.5 

-23.4 

-1.2 

-0.3 

-0.9 

-220.7 
10.1 
2.1 
0.8 
3.3 
0.0 
-298 



as 



contribution of the remainder, while small, is not negli- 
gible for n < 20 and is calculated in the random-phase 
approximation (RPA) for these states. It is essentially 
zero for larger n and is evaluated in the DHF approxi- 
mation for n > 20 without loss of accuracy. We find that 
it is essential to use RPA approximation to evaluate the 
tail contribution for as as DHF calculation significantly 
overestimates the tail yielding —25 Og while RPA calcu- 
lation gives —16 Oq. DHF approximation is expected to 
significantly overestimate the value of the tail contribu- 
tion as it significantly overestimates the value of the main 
term. The autoionizing contribution is very small, 0.1%, 
for as and is negligible at the present level of accuracy 
for j3s- It is evaluated in the DHF approximation. 

In Tables |2 and ^I we present the details of our as 
and Ps calculations for the transition between the ground 
state 6P1/2 and the excited 7Pi/2 state in Tl I. We sepa- 
rate the contributions from the nS and nD^/2 terms given 
by the Eqs. I^I^J. Furthermore, the contribution from 
each n in the Eq. |5j for the main term is given sep- 
arately in the corresponding (nlj) row to demonstrate 
rapid convergence of the sums in Eqs. Q] ^. The tail 
contributions for nS and nD^/2 sums are listed in the 
rows following the corresponding main term. The total 
contribution from the core-excited autoionizing states are 
listed in rows labeled (5a|°" and (5/3|°'<= of Tables U and 
im respectively. 

We list the values of the reduced electric-dipole ma- 
trix elements and energies used in the calculation of as 
and Ps in Tabled Since the same matrix elements and 
energies contribute to as and f3s, we do not repeat the 
energy and matrix element values in Table |n] We use 
recommended values of energies from the National Insti- 
tute of Standards and Technology (NIST) database [ij 
when calculating all of the polarizability values in this 
work. The corresponding energy differences are listed in 
columns two and three of Table ^ Electric-dipole ma- 
trix elements evaluated using the SD all-order method 
(Eq. 0) are given in columns labeled Z^^. It should be 



Ws 



This work 

Expt." 

Expt.'' 



368 
247±12 

377±8 



298 

198±10 

313±8 



noted that the values of Z^^^ given in Table ^ generally 
differ by 0.4-1.0% from the values of Z^^ presented in 
Ref. 0. These differences arise because we include the 
Breit interaction here and, to the lesser extent, because 
of more accurate basis set used in the present work. 

We also conducted a semiempirical scaling procedure, 
described for example, in Refs. |3,|i3|, for three transi- 
tions, 7S - 7Pi/2, 7S - 7P3/2, and 7Pi/2 - 6D3/2. These 
transitions give significant contributions to the as and 
(3s as well as Stark shifts that we discuss below. The 
scaling is carried out by multiplying the values of the 
corresponding single valence excitation coefficients by the 
ratio of the theoretical and experimental correlation en- 
ergies and repeating the calculation of the matrix ele- 
ments with the modified excitation coefficients. Study 
of the breakdown of the correlation corrections showed 
that the dominant contribution to the correlation correc- 
tion to the values of these matrix elements comes from 
a single term that contains only single valence excitation 
coefficients. It has been shown that the scaling works 
effectively in such cases (see, for example [i3il3) since 
it is specifically aimed at correcting the dominant con- 
tribution. In the case of the 75* — 7P transitions, the 
correlation breakdown is the same as in Cs, where the 
scaled values are in excellent agreement with the high- 
precision experiment. The scaling can not be applied 
to improve the values of the 61)3/2 — 6P1 /2 matrix ele- 
ment as the the corresponding correlation correction is 
not dominant, leading to possibly reduced accuracy of 
this transition in comparison to 75" — 7Pi/2, 7S — 7P3/2, 
and 7Pi /2 — 61)3/2 ones. 

We find that the contribution of the 6D3/2 term from 
Eqs. 1^ 121) dominates the values of as and /3s. While 
all valence terms contributing to as have the same sign, 
there is significant cancellation of terms contributing to 
(3s- In fact, all nS terms contribute with the sign op- 
posite to that of the ^1)3/2 terms with the exception of 
the first ones, 75" and 6D3/2. Owing to this cancellation, 
our value of Ps niay be somewhat less accurate than the 
value of as- 

In Table IIIII we compare our results for the Stark- 
induced scalar and vector polarizabilities for the 6P1 /2 — 
7Pi/2 transition in Tl with experimental measurements 
by DeMille et al. [j] and Tanner and Commins Ql- The 
conversion factor between the units of (hb/c) (cm/V) 
used in Refs. [l|, Q and atomic units used in the present 
work is 10~^a£;/,/2eao=l. 8762x10^, where Eh is Hartree 



energy. Our results support the measurements of DeMille 
et al. [l| and clearly disagree with the measurements re- 
ported by Tanner and Commins [3| ■ Our value of as is 
nearly within the experimental uncertainty of the Ref. 'l] 
measurement and our value of Ps differs from the central 
experimental value of Ref. Ij by 2a. As we discussed 
above, accurate calculation of (3s is more difficult owing 
to the cancellation of the nS and nD^/2 contributions. 



III. STARK SHIFT WITHIN THE 6P1/2 - 7Si/2 

AND 6P1/2 - 7Pi/2 TRANSITIONS IN ATOMIC 

THALLIUM 



We calculate the Stark shifts within the QP1/2 — 751/2 
and 6P1/2 — 7Pi/2 transitions as a differences of scalar 



dipole polarizabilities a of the 6P1/2 ground state and 
the 75*1/2 or lPi/2 excited states. The expression for a 
is given by (see, for example, Refs. |5L[l5|): 



TABLE IV: The contributions to the scalar dipole 6P1/2 and 
7P1/2 polarizabilities Qf(a.u.) in Tl. 



Contribution 


nlj 


q(6Pi/2) 


«(7Pl/2) 


a'^''''\nS) 


75i/2 


9.2 


-345 




8<S'i/2 


0.5 


616 




9Si/2 


0.2 


13 




105i/2 


0.1 


3 


a*^"(n5') 




0.7 


6 


a-^'"(nD3/2) 


6-D3/2 


11.0 


4280 




7D3/2 


2.1 


217 




8£>3/2 


0.7 


39 




9P>3/2 


0.4 


14 


a*^"(n^3/2) 




4.7 


47 


^corc 




24.1 


24 


(5q^°'° 




-3.3 





Total 




50.4 


4915 



a{nQPi/2) = ^[ls{nSi/2) + Is{nD3/2)] 

n 

a{noSi/2) = ^[ls{nPi/2)+Is{nP3/2)] : 



where 



Is (nlj) = 



1 



y2 

nol.ojo,nlj 



3 Enlj — Enol, 



ojo 



Our results for the 6P 1/2, IP1/2, and 75i/2 polarizabili- 
ties are given in Tables iTVl and where we use the same 
designations as in the previous section. The polarizabil- 
ity Qicorc of the Hg-like ionic core is also evaluated using 
RPA approximation. A more detailed discussion for the 
encore in Na, K, Rb, Cs, and Fr atomic systems is found 
in Ref. 7] . We note that contributions from Ofcorc cancel 
when we evaluate the Stark shift for a transition. 

Only the polarizability contributions are listed in Ta- 
ble ^3 since all relevant energies and electric-dipole ma- 
trix elements are already listed in Table Q] The calcula- 
tion of the 75i/2 polarizabilities involves the calculation 
of other series of the matrix elements, 75 — nPi/2 and 
75 — nP-ii2- We Hst those values, calculated using SD 
all-order method, together with the corresponding en- 
ergy differences taken from [ij in Table IVl 

The value of the 6P1/2 polarizability has two domi- 
nant (and nearly equal) contributions, from6Pi/2— 6Z?3/2 
and 6P1/2 — 75 transitions. The value of the 7Pi/2 
polarizability is dominated by the contribution from 
7Pi/2 — 6133/2 transition. The value of the 75 polarizabil- 
ity is dominated entirely by the contributions from both 
IS — IP transitions. As we have discussed above, we con- 
ducted more accurate calculation of the 7Pi/2 — 6I33/2, 
75— 7P1/2, and 75— 7P3/2 electric-dipole matrix elements 
by rescaling the single valence excitation coefficients with 
the correct value of the correlation energy leading to more 
accurate evaluation of the dominant contributions to 75 



TABLE V: The contributions to the scalar dipole IS polariz- 
ability a in Tl. The energy differences and absolute values of 
the electric-dipole reduced matrix elements for relevant tran- 
sitions are also listed. All values are given in a.u. 



Contribution 


nlj 


Enlj — Ets 


\Z7S,nlj\ 


«(7Si/2) 


a--°(nPi/2) 


6P1/2 


-0.120640 


1.826 


-9.2 




7Pl/2 


0.035004 


6.016 


344.7 




8P1/2 


0.067847 


0.706 


2.5 




9Pl/2 


0.081574 


0.296 


0.4 


«*""(" A/2) 








0.4 


a-^-(nP3/2) 


6P3/2 


-0.085134 


3.397 


-45.2 




7P3/2 


0.039565 


8.063 


547.7 




8P3/2 


0.069545 


1.474 


10.4 




9P3/2 


0.082401 


0.713 


2.1 


«*""(" A/2) 








3.0 


^corc 








24.1 


Total 








880.8 



and 7Pi/2 polarizabilities. Our value for the a(6Pi/2) = 
50.4 ao'^ is in good agreement with the theoretical result 
49.2 ao^ by Kozlov et al. [H]. 

In Table IvTI we compare our results for the Stark shift 
in the 6P1/2 — 75i/2 and 6P1/2 — IPin transitions with 



TABLE VL The all-order values of the ground and excited 
state polarizability differences are compared with experimen- 
tal results by Doret et al. Q - a, by DeMille et al. Q- b, and 
by Fowler and Yellin iq- c. 





a(6Pi/2) - 


a(7Si/; 





q(6Pi/2) - 


«(7Pi/2) 


Present work 


-830 






-4866 




Expt." 


-829.7±3.1 










Expt.'' 


-900±48 






-4967±249 




Expt.'^ 


-776±80 











experimental results from Refs. 0,yil3- The conversion 
factor between the lS.vs in kHz/(kV/cm)^ units used in 
Refs. jH, 01 0| to polarizabilities in atomic units used in 
the present work is 2 x 10~^/i/(47reoao) = 8.03756, where 
h is the Plank constant. Our result for the 6P1/2 — 75i/2 
Stark shift agrees with the most accurate, 0.4%, exper- 
imental value from Ref. Q within the experimental un- 
certainty. The all-order value of the Stark shift for the 
6^1/2 ~ TPi/2 transition also agrees with the experiment 
within the experimental uncertainty. We note that this 
comparison essentially tests the accuracy of the 75' and 
7Pi/2 polarizability calculations since the ground state 
polarizability is small comparing to the excited state po- 
larizabilities. 



75i/2 and &P1/2 ~ TPi/2 transitions and the Stark- 
induced amplitudes as and /3s in the 6P1/2 ~ 7Pi/2 tran- 
sition in atomic thallium is presented. Our results a = 
368 ao"^ and |/3s|= 298 ao"^ support the experimental mea- 
surements carried out by D. DeMille, D. Budker, and E. 
D. Commins [Phys. Rev. A50, 4657 (1994)]. Our results 
for Stark shifts in the 6P1/2 — 7iS'i/2 and 6P1/2 — 7Pi/2 
transitions agree with the most accurate measured val- 
ues within the experimental uncertainties for both tran- 
sitions. 
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